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New coherent states and a new proof of the Scott correction 

Jan Philip Solovej and Wolfgang L Spitzer 

Abstract. We introduce new coherent states and use them to prove semi- 
classical estimates for Schrodinger operators with regular potentials. This can 
be further applied to the Thomas-Fermi potential yielding a new proof of the 
Scott correction for molecules. 



1. Introduction 

In this paper we review a novel proof of the Scott correction for neutral mol- 
ecules. So suppose, we have M nuclei of positive charges Z = {Zi, . . . , Zm) G IR.+^ 
located at positions R — {Ri, . . . ,Rm) G M.^^ ■ We choose the charge of an electron 
equal to —1, so that neutrality is expressed as \Z\ = J2j=i — where N is the 
number of electrons. Further, we use atomic units where It? = m. 

The interaction of a single electron with all the nuclei is equal to 

M ^ 

(1.1) n^^^--) = E^- 

We now write the molecular non-relativistic Schrodinger operator in the form 
H{Z, R) — H{Zi, . . . , Zm',Ri, • • ■ , Rm) 

iX^ViZ,R,x,))+ J2 T^^r 

The operator H{Z, R) acts on the space Aili L'^i^^ x {—1, 1}), where ±1 refers to 
the spin variables. We are interested in the ground state energy, 

(1.2) E{Z,R) ^MspecH{Z,R), 

and in particular, in an asymptotic expansion for large charges. Let us state the 
main theorem in this paper. 
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Theorem 1.1 (Scott correction). Let Z = \Z\{zi, . . . ,zm), where zi,...,zm 
> and R = \Z\~^/^{ri, . . . , tm), with \ri — rj \ > tq > 0, for all i ^ j. Then, 

(1.3) E{Z,R) = E^^{Z,R) + ^ +0(|^|'"'/'°), 

l<j<M 

as \Z\ — > oo, where the error term 0{\Z\'^~^^^^) besides \Z\ depends only on 
zi,...,zm, andro- 

The leading Thomas-Fermi (TF) term, which is of the order IZ]"^^^ was first 
rigorously derived in the seminal work by Lieb and Simon [LSi] using the Dirichlet- 
Neumann bracketing method. 

The Scott correction, i.e., the term ^ 5^i<j<M -^j '^^'''^^ proven by Hughes [H] 
(a lower bound), and by Siedentop and Weikard [SW] (both bounds) in the case of 
atoms. The atomic case is simpler since in TF theory atoms are spherically sym- 
metric. Bach [B] proved the Scott correction for ions. Finally, Ivrii and Sigal [IS] 
accomplished a proof of the Scott correction for molecules, which was recently ex- 
tended to matter by Balodis Matesanz [M]. Here, we present another proof for 
molecules. 

It was later shown by Lieb [LI] (and independently by Thirring [T]) how co- 
herent states can be used to give a simple proof of the leading TF term with good 
upper and lower bounds; sec also a recent improvement by Balodis Matesanz and 
Solovej [MS] . We want to stress that in order to prove an asymptotic expansion for 
E{Z, R) capturing the Scott term one basically needs to prove a local trace formula 
for regular potentials (see Theorem 4.1 with n = 3) up to the order where e 

is any positive number. We accomplish e = 1/5. 

A quick explanation for the Z^-correction goes as follows. Whereas the leading 
TF term comes from the bulk of electrons, the correction comes only from electrons 
close to the nuclei where the Coulomb attraction is unscreened by the presence of 
the other electrons. From the exact solution of the hydrogen atom one may extract 
the Scott correction (sec [LI]). Notice, that the Scott correction for molecules is 
just the sum of the corresponding atomic corrections. This is not the case for the 
leading term. 

This review is organized as follows. In Section 2 we recall the main analytic 
tools and state the main properties of the TF potential. We introduce the new 
coherent states in Section 3. In Section 4 we sketch the proof of the semi-classical 
estimates on the sum of negative eigenvalues for regular and the TF potential. In 
the last Section we present the proof of the main Theorem 1.1. For more details 
we refer to our paper [SS]. 

2. Preliminaries 

2.1. Some Inequalities. Here we collect the main inequalities which we need 
in this paper. Various constants are typically denoted by the same letter C, and in 
all cases sharp constants do not play a role. 

Let p> 1, then a complex- valued function / (and only those will be considered 

here) is said to be in Lp(IR") if the norm := (/ |/(.t)|p rf.x) is finite. For 

any 1 < p < t < q < oo we have the inclusion D L'^ C L', since by Holder's 
inequality < ||/||^||/|| J"^ with Xp-^ + (1 - X)q-^ = t-\ 
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We call 7 a density matrix on i^(]R") if it is a trace class operator on iy^(M") 
satisfying the operator inequality < 7 < 1. The density of a density matrix 7 is 
the function p-y such that Tr(7^) = / p^{x)9{x)dx for all 9 e Co°(K") considered 
as a multiplication operator. 

If € L^(M^ X {—1, 1}) is an A''-body wave-function, then its one-particle 

density, p^, is defined by 

N 

Pi,{x)=^ X! "' X] i \i^{xi,si;...;xN,SN)\^S{xi-x)dxi---XN- 

i—l si— ±1 s^v— il 

The next inequality we recall is crucial to most of our estimates. 

Theorem 2.1 (Lieb-Thirring inequality) . One-body case: Let 'y be a 

density operator on L^(R"), then we have the Lieb-Thirring (LT) inequal- 
ity 



(2.1) TV[-iA7] >Jr„yp, 



1+2/n 
7 ' 



where Kn is some positive constant. Equivalently, let V € i^+"/^(R") 
and 7 a density operator, then 

(2.2) Tr[(-iA + l/)7] > -L^ j \V_\'+^'\ 

where X- := min{x, 0}, and Ln some positive constant. 
Many-body case: Let ip e A^i L'^i^^ x {-1, !})• Then, 

(2.3) ^^,f:-iA,V')> 2-2/^X3 Jpf. 



The original proofs of these inequalities can be found in [LT]. From the min- 
max principle it is clear that the right hand side of (2.2) is in fact a lower bound 
on the sum of the negative eigenvalues of the operator —^A-\-V. 

We shall use the following standard notation for the Coulomb energy: 

D{f) = D{f,f) = 1 ff f{x)\x - y\-'f{y)dxdy. 



2 . . 

It is not difficult to see (by Fourier transformation) that ||/|| :— D{f Y^^ is a norm. 

Theorem 2.2 (Hardy-Littlewood-Sobolev inequality). There exists a constant 
C such that 

(2.4) D{f)<C\\f\\l„. 

The sharp constant C has been found by Lieb [L4], see also [LL]. 

Finally, we state the two inequalities which we shall need to estimate the many- 
body ground state energy, E{Z, R), by an energy of an effective one-particle quan- 
tum system. The first one is the electrostatic inequality providing us with a lower 
bound. This inequality is due to Lieb [L3], and was improved in [LO]. 

Theorem 2.3 (Lieb-Oxford inequality). Let ip € L^(K^^) be normalized, and 
Pif, its one- electron density. Then, 

-1,/, \ Ti(^.\_n I 4/3 




(2.5) (V, 1, |x,-x,r>) >I)(p^)-C / 
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An upper bound to E(Z, R) is furnished by a variational principle for fermionic 
systems. This is also due to Lieb [L2]. 

Theorem 2.4 (Lieb's Variational Principle). Let j be a density operator on 
L^(R^) satisfying 2Tr7 = 2 J pj{x) dx < Z (i.e., less than or equal to the number 
of electrons) with the kernel p-yix) = ^{x, x). Then, 

(2.6) E{Z,R) < 2TT\{-^A-ViZ,R,x))j\ + D{2p^). 

The factors 2 above are due to the spin degeneracy. The ground state wave- 
function carries a spin and is really a function on L^(1R.^^; ), but only its spatial 
dependency is of interest here. 

2.2. Thomas-Fermi Theory. Here we quickly state the properties about TF 
theory which are needed for our proof. The original proofs can be found in [LSi] 
and [LI]. 

Consider z = (2:1, . . . , 0m) e and r = (n, . . . ,rM) e K^*^. Let < p e 
L^/^{M?) n L1(IR3) then the TF energy functional, S'^^ (omitting the nuclei-nuclei 
interactions), is defined as 

(2.7) f TP(p) = ^(377^)2/3 j p{xfl^ dx-j F(z, r, x)p{x) dx + D{p), 

where V is as in (1.1). 

Theorem 2.5 (Thomas-Fermi minimizer). For all 7, = {z\, . . . , zm) € and 

r = (ri, . . . jTm) G IR"^^ there exists a unique non-negative p'^^{z,r,x) such that 
J p^^{z, r, a;) dx = Y^^L-^ zu and 

fTF(pTF^^i^f|^TF(^) : 0<pei5/3(^3)^j^l(^3)|_ 

We shall denote by E'^^{z,r) := £'^^{p'^^) the TF-energy. Moreover, let 

(2.8) V^^iz, r, x) := V{z, r, x) - p"^^ * |a;|-^ 

be the TF-potential, then V'^^ > and p^^ > 0, and p^^ is the unique solution in 

L^/'i{U:i)f}L^{R^) to the TF-equation: 

(2.9) yTF(z, r, x) - i(37r2)2/3pTF(^^ ^ 



2/3 



Very crucial for a semi-classical approach is the scaling behavior of the TF- 
potential. It says that for any positive parameter h 

(2.10) V'^^{z,T,x) = h-^V^^ih\h-^r,h-^x), 

(2.11) pT^(z,r,.x) = h-'^p'^^{h^z,h-h,h~^x), 

(2.12) ^^^(z,r) = h'^E^''{h^z,h-h). 

By h~^r (and likewise for h^z) we mean that each coordinate is scaled by h~^. 
Notice that the Coulomb-potential, V, has the claimed scaling behavior. The rest 

follows from the uniqueness of the solution of the TF-energy functional. 

We shall now state the crucial estimates that we need about the TF potential. 

Let 

(2.13) d{x) =min{|a;-rfe| | k = 1,...,M}, 
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and 

(2.14) f{x) = inin{d{x)-^/^,d{x)-^}. 
For each k = 1, . . . ,M we define the function 

(2.15) Wk{z,r,x) = V^^{z,r,x) - Zk\x-rk\-^. 

The function Wk can be continuously extended to a; = rfe. We have the following 

estimate for the TP potential. 

Theorem 2.6 (Estimate on TP potential). Let z = {zi,...,zm) e ^+ and 
r = (ri, . . . , tm) € M.^^ . For all multi-indices a and all x with d{x) ^ we have 

(2.16) \d^V^^iz,r,x)\<C^f{xfd{x)-\''K 

where Ca > is a constant which depends on a, zi,. . . , zm, and M. 

Moreover, for \x — rk\ < rinin/2, where rmin = min^^^ \rk — re\ we have 

(2.17) 0<Wk{z,r,x)<Cr-l + C, 
where the constant C > here depends on zi,. . . , zm, and M. 

The relation of TF theory to semi-classical analysis is that the semi-classical 
density of a gas of non-interacting electrons moving in the TF potential is simply 
the TP density. More precisely, the semi-classical approximation to the density 
of the projection onto the eigenspace corresponding to negative eigenvalues of the 
Hamiltonian — — V'^^ is 

2 / l-^=2'/'{3n')-'{V^^r/^z,r,x)=p^^{z,r,x). 

Here, the factor two on the very left is due to the spin degeneracy. Similarly, the 
semi-classical approximation to the energy of the gas, i.e., to the sum of the negative 
eigenvalues of — |A — V'^^ is 

(2.18) 2/(ip^-y-(z,r,.))_^ = ^"'(-'r,.)V^.^. 

= E'^^iz,r)+D{p'^^). 
In Section 4 we shall make the semi-classical approximation more precise. 

3. New coherent states 

Coherent states provide a natural semi-classical description of quantum me- 
chanics. We shall denote these states by \u,q), where {u,q) G M^" is a point in 
phase-space. Their wave-function is given by 
(3.1) {x\u, q) = (7r/i)-"/*e-(^-"''/2/»gWft_ 

Let Tiu,q = \u,q){u,q\ be the projection onto the coherent state \u,q), then they 
satisfy the completeness condition (in the sense of quadratic forms) 



(3.2) J n„,, 



dudq ^ 



''"^ (27r/i.) 

As functions on phase-space the coherent states are localized on a scale of the order 
of h. We want to broaden this by defining the operator 

(3.3) Qu,q ■= j w{u-u',q- q')Uu',q' du'dq' 



6 



JAN PHILIP SOLOVEJ AND WOLFGANG L SPITZER 



with 



\7r(l - ha) J 



The new scale is 1/a > h, which becomes clearer when we look at its kernel, 

(3.4) gu,g{x, y) = (7r/i)-"/2e-''('^-'')'+^«(^-^)/^-iS^(^-^)'. 

For simplicity, we have chosen a Gaussian weight, w, in the definition of Qu,q- 
We shall use the operators Qu,q as our new coherent states. ^ Note that if we 
let a — > l//i then Qu,q converges to n„_g. A straightforward calculation gives the 
following result. 

Lemma 3.1 (Completeness of new coherent states). These new coherent oper- 
ators satisfy 

This resolution of the identity provides us with a representation of Schrodinger 
operators as phase-space integrals. This will be useful when we prove a lower bound 
on the sum of the negative eigenvalues of Schrodinger operators. 

Theorem 3.2 (Coherent states representation). Consider functions F and V 
in C^(R"), for which all second and third derivatives are bounded. Let a{u,q) = 
F{q) + V{u), then we have for a < 1/h and b = 2a/(l + /i^o^) the representation 

(3.6) F{-ihV)+V{x) = J Qu,qHu,qgu,q-^^+^ 
as quadratic forms on Cq°(R") vdth the operator-valued symbol 

(3.7) Hu,q = a{u, q) + ^Ac7(u, q) + 9„cr(u, q){x - u) -\- dq<j{u, q){-ihV - q). 
The error term, E, is a bounded operator with operator norm 

(3.8) ||E|| < l|3V||oo + Ch% J2 ll^"^lloo- 

|a|=3 |a|=2 

Starting with the identity (3.5), the representation of Schrodinger operators as 
in (3.6) arises by splitting the product ^ apart while sandwiching the symbol 

Hu,q- This operator- valued symbol can be thought of as the first order Taylor 
expansion of the classical symbol cr{u, q) at (x, —ihV). Clearly, one could consider 
higher order expansions but this in not needed here. Also notice that as a | 1/h the 
linear term in (3.7) does not contribute in (3.6) and one gets the familiar classical 
approximation a{u, q) + i^IS.a[u, q). The representation (3.6) is symmetric in space 
and momentum due to the symmetric Gaussian weights in the definition of 



Sometimes (e.g. see [LI]) it is useful to consider other coherent states where the Gaussian 
function in (3.1) is replaced by a general function. Similarly, one could use them to define 
corresponding generalized coherent states but from a computational point of view the above choice 
is the simplest. 

These new coherent states should not be confused with the quantum coherent operators 
introduced by Lieb and Solovej in [LSo] in order to compare two quantum systems. 
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One major advantage of coherent states is that a positive (upper) symbol leads 
to a positive operator. This is important when writing down explicit variational 
states and brings us to consider more generally operators of the form 

(3.9) j gu,qf{Au,q)Qu,gdudq. 

Here, Au,q — Bq{u, q) + Bi{u, q) ■ x — ihB2{u, (jf) • V is a Hcrmitian operator which 
is linear in x and —ih\7^ and / : K ^ M is any polynomially bounded real function. 
We shall denote by A^.q the linear function Au,q{v,p) = _Bo(w. <?) + Bi{u.q) ■ v + 
B2(u,q) - p. When Au.q{v,p) is independent of (w,p), i.e., if Bi = B2 — and if 
a ^ we recover the usual coherent states representation of an operator. Thus 
on the one hand, we do not use as sharp a phase-space localization as the one- 
dimensional coherent state projection since a < 1/h, but on the other hand, we use 
a better approximation than if Au^q were just a constant. 



4. Proof of semi-classical estimates 



4.1. Regular potentials. The key application of coherent states will be a 
proof of a semi-classical expansion of the sum of negative eigenvalues of (local- 
ized) Schrodingcr operators. Wc shall restrict ourselves to localization fimctions 
supported in balls. Recall that wc use the convention, X- = (a;)_ = min{.T,0}. 

Theorem 4.1 (Local semi-classics). Let n> 3, (j) G Cq~^'^{W^) be supported in 
a ball El of radius £ > 0. Let V e C^{Bi) be a real potential. Let H = -/i^A -|- V, 
h> and a{u, q) = q^ + V{u). Then for all h> and f > we have 



(4.1) 



Ti[<j>H(P]^ - (27r/i)" 



(u)^cr(w, q)- dudq 



where the constant C depends only on the dimension n, 

(4.2) sup ||£l"l5"<^||oo, and sup Wf-^l^'^^d'^VW^. 

|a|<n-|-4 |a|<3 

Moreover, there exists a density matrix 7 such that 

(4.3) Tr[(/.iJ(/.7] < (27r/i)-" j (l){ufa{u, g)_ dudq + C/j-"+6/5_^n-F4/5^n-6/5^ 
and such that its density Pj{x) satisfies 

(4.4) p^{x) - (27r/i)-"a;„ \V{x)-f/^ < Ch-^+^'^'' f^-^'^H-^'^'' , 



for (almost) all x G B(, and 
(4.5) 



j (l){xfp^{x)dx-{2'Kh)-''ujnj (t>{xf\V{x)-\''''^ dx 



^ — n-t-6/5yn— 6/5^n— 6/5 



The constants C > in the above estimates again depend only on the dimension n, 
the parameters in (4-2), and the volume of the unit ball in M", ujn. 

As mentioned in the Introduction, any power 0(/i^"+^+^) with e > is suffi- 
cient to prove the Scott correction in the main theorem (1.1). The power | comes 
from optimizing the error bound in (3.8) by choosing 6 = 
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Sketch of proof. By a simple scaling argument we may restrict ourselves 
to the unit ball setting £ = 1 and the case / = 1. We start with a sketch of the 
lower bound. We may assume that h is sufficiently small. Using the representation 
(3.6) we have that 



Tr[(j)H(t)\- > Tr 

(4.6) +Tr 

Here, < e < 1/2, and 

Hu,g = ^{u, q) + ^Act(u, q) + dua{u, q){x - u) + dqa{u, q){-ihW - q) 

with a{u, q) = [1 — e)q^ + V{u). Utihzing the variational principle for the first trace 
and the LT inequality for the second one we obtain the bound 



(27r/i)"Tr[0i70]_ > / Tr 



7u,q 



H 



u.q 



dudq - C£-"/2(6-3/2 + hHy+^/\ 



We shall eventually choose e = j{b^'^^^ + h^b). Since Hu^q is a linear combination of 
X and V this operator can be easily rotated into the momentum operator alone (up 
to some constant term). Then, we conveniently have an expression for the negative 
part of Hu^q, and a fortiori, the trace becomes a Gaussian-like integral which we 
have to estimate. In this integral, the linear function 

Hu,q{v,p) = a{u,q) + ^Aa{u,q) + dua{u,q){v - u) + dqa{u,q){jp - q) 

replaces the operator kernel of -ffu.q; notice that (x, — i/iV) is simply substituted by 

We can show that if we consider the u-integration over u outside the ball of 
radius 2 then this is bounded below by —Cb~^/'^. On the other hand, the integration 
over B2 can be estimated from below by 

/ (j)^ {y + h^ab{u — v)) Gb{u — v)Gb{q — p) [Hu,qiv,p)]_ dudqdpdv, 

with Gb{v) = (5/7r)"/2exp(-5w2). Since we are looking for a lower bound we may 
as well extend the last integral to K". Notice that we may now perform the p- 
integration and obtain 



(27r/i)"Tr[<?!)i?0]- 



> 



'n + 2 



(!-£)- 



0^ (v + h^abu) Gb{u)Gb{q) 



V{v)+Uu,q)-C\u\{b-' + \u\^) 
(4.7) -C(6-3/2 ^ f^2^-^^ 

where we have introduced the function 

iv{u,q) = ^Aa{v,Q) - (1 - e)g^ - ^^didjV{v)uiUj 



r + l 



dudqdv 



HI 
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By expanding wc find that 

V{v)+l{u,q)-C\u\{b-' + \u\^) 
< \V{v)\I+' + Q + l)\V{vUnviu,q) 

+c(\Uu,q)\+C\u\{b-^ + \u\^y 
+C|«|(6-i + 

We have used that since n > 3, the function x i— > la;-!^"*"^ is C^(M"). Hence, 

2w„ 



(27r/i)"Tr[(/.if0]_ 



n + 2 



{l-sY 



02 („ + h'^abu) Gb(u)Gbiq) 

+ + l) \V{v)-\^Uu,q)) dudqdv 
-C{b-^'^ + h'^b). 

We now expand (jp' and use the crucial identities for Gaussian integrals, 
j ^v{u,q)Gb{u)Gb{q)dudq = and J uGb{u)du = 0. 



We arrive at the lower bound, 

2w„ 



{2nh)"-Ti[(j)H(f)l 



> — 



n + 2 



(l-£)-^ J (j){vf\V{v)\I^^dv - C{b-^/^ + hH) 
{vfa{v,p)-dvdp - C{b-^'^ + h%). 



Finally, we choose a = h^^^^ and e = b^^^^. 

Now we come to the upper bound. We shall show here only the construction 
of the density matrix 7. Let x = X(-cx3,o] be the characteristic function of (—00, 0] 
and 



_ ( a{u, q) + ^Acr(w, q) + 9„cr(u, q){x - u) + dga{u, q){-ihSI - q) if u e B2 
~\ iiu^B2 

We then define 

(4-8) 1 = jQu,,x[hu,,]Qu„^. 

Since < x[^«,9] < 1 it is obvious that < 7 < 1. The arguments showing that 

Tr[7(/)i?0] < {2-Kh)-'' j 4)'^{u)a{u,q)-dudq + Ch-''+^l^ 
are then very similar to the above calculations for the lower bound, see [SS]. □ 

4.2. Thomas-Fermi potential. In this Section we shall sketch the proof of 
the Scott correction for the TF potential. 
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Theorem 4.2 (Scott corrected semi-classics). For all h> and all n, . . . ,rM 
e M.^ with mirife^TO \rm — rk\ > ro > we have 

M 

Tr[-/i2A-l/TF]__(27r/i)-3 



(4.9) 



fe=i 



where C > depends only on z\,...,Zm, M, and tq. Moreover, we can find a 

density matrix 7 such that 

(4.10) Tr [{-h^A - yTF)^] < Tr [~h^A - V^^] _ + C/i-^+Vio, 
and such that 

(4.11) D(^p,-^iv--r/^)<ch-^^^/\ 

and 

(4.12) y < ^ I V^^{xf/'dx + Ch-'+'/\ 

with C depending on the same parameters as before. 

Sketch of proof. From Theorem 2.6 we know that the TF potential has an 
inverse fourth power law decay far from the nuclei. Thus, a region outside some 
ball of radius R (which scales with h) should contribute little to the sum of negative 
energies. For this purpose we introduce a first partition of unity. So let us choose 

(4.13) i? = /i-i/2, 
and consider functions $j- S C°°(]R") such that 

(1) = 1, 

(2) $_(a;) = 1 if d(x) < R and $_(x) = if d{x) > 2R. 

Denote T = (V$_)^ + (V$+)^. Then T is supported on a set whose volume is 
bounded by CR^ (where as before C depends on M) and ||I||oo < CR~^. Using 
the standard IMS localization formula and then the LT inequality we find that 

> TV[$_(-/i2A - _ /j2j)$_j_ + Tr[$+(-/i2A _ yTF _ h''l)<^+\_ 

> TV[$_(-/i2A - V^^ - /i'j)$_]_ - Cih'R-^ + h-^R-'). 

With the chosen R the last term is of the order /i^^/^. 

On the support of we want to use the hydrogenic approximation of the TF 
potential close (of the order of h) to the nuclei, and on the rest the semi- classical 
estimates from the previous Section. Let us introduce the function 

M _^ 

(4.14) £(^)= l(^l+^(|a;_r,|2+/j2)-l/2j . 

fc=l 

Note that f is a smooth function with 

0<e{x)<l, and ||V^(x)||oo < 1. 

Now, we fix some localization function (j) G C|^(M'^) with support in the unit ball 
{|a;| < 1} and such that / (j){xfdx = 1. It is not difiicult (cf Theorem 22, [SS]) 
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to find a corresponding family of functions (j)u € C^(M'*), u e M'', where (pu is 
supported in the ball {\x — u\ < £{u)}, with the properties that 



(4.15) J Mxf£{u)-^du = l and ||5"<^„||oo < C^)" 



|a| 



for all multi-indices a, where C > depends only on a and </>. 

One can show from (2.16) in Theorem 2.6 that for all u e M" with d{u) > 2h, 



(4.16) 



sup \d"v^''{x)\<cf{ufe{u)-\'^\, 

\x-u\<e(u) 



where C > depends only on a, zi, . . . , zm, and M . This is the requirement for 
the semi-classical estimates from Theorem 4.1 to apply with i{u) i, f{u) f. 
Another application of the IMS formula shows that 

(4.17) Tr[-/i2A - 

> / Tr[<^„ (-/i^A - V^^ - Ch''£{u)-^) ^u]-i{u)-^du - Ch-^'^. 

Jd{u)<2R+l 

By similar arguments we get corresponding estimates for the hydrogenic operators 



replacing V by 



\x-rk\ 



1 in the above estimates. In particular, if we choose h so 



small that R > max^l^fc} then on the support of we have —Zk\x — rk\ ^ + 1 > 0. 
Thus we have 



(4.18) Tr 



Zk 



> 



\x-rk\ 
Tr 



+ 1 



Jd(u)<2R+l L ^ 



ld{u)<2R+l 



Zk 



\x-rk\ 



+ 1 



We arrive at analoguous upper bounds if we utilize the density matrix 

7 = / cP^x{M-h^^-V^^)^u)Mu)-^du 

Jd{u)<2R+l 

as a trial operator. I.e., 

(4.19) = J Tr[(l)u{-h''A-V^'')(l)u]-e{u)-^du. 



d{u)<2R+l 



Similarly, 
(4.20) Tr 



Zk 



+ 1 



x-rk\ 
< I Tr 

d{u)<2R+l 



-fi'A 



Zk 



\x - rk\ 



i)<Pu\j{u)- 



'du. 
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We now introduce the quantities 



M 



2fe 



fc=i 



\x-rk\ 



+ 1 



M 

£»_(«) := ^TrU„(-/i2A-_i^ + l-C7/i2^(u)-2),^„ 



and 



1 

Dsc(w) := {2nh)-^ J Mxfip''-V^^ix))-dpdx 



M 



\X - Tfe 



+ 1 



Then, from (4.17), and (4.19) we have 



(4.21) 



M 

TV[-/i2A - yTF]_ _ V Tr [-/i^A - + 1 



fc=i 



> 



/ n+(^^)£(^^)-3du-C/l-^/^ 



and from (4.15) we get 

/M „ 



Zk 



1 ) (ipda; 



/s=l" 



(4.22) 



\x-rk\ 
= j Dsc{u)eiu)-^du. 



Next, we compute explicitly both the quantum and the semi-classical energies for 
the Coulomb potential, namely 



Tr 
and 



Zk 



\x-rk\ 



( 4/j2 



{2nh)-^J {p- 



l<n<Zk/{2h) 
Zk 



+ 1 



^ dudp ■■ 



\u-rk\ ' V-"""^ 12/i3" 

The first statement of the theorem is thus proven once we establish lower bounds 

on D+{u) — Dsc{u) and D^{u) + Z?sc(w). Here, we have to distinguish between 
the region d{u) < 2h and the semi-classical region, 2h < d{u) < 2i? + 1. 

In the region close to the nuclei, d{u) < 2h, we use the estimate (2.17) on the 
potential, Wfe(x) ~ V'^^{x) — Zk\x — rk\~^- Then, all bounds on £>± (w) are obtained 
via the LT inequality. 

Secondly, let 2h < d{u) < 2R + 1. On the ball {x \ \x - u\ < ({u)}, the TF 
potential satisfies the estimate (4.16) and satisfies (4.15). Hence, we may use 
Theorem 4.1. A similar semi-classical estimate holds for the Coulomb potential, 

\x-rk\ 
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The density matrix which we choose for (4.10-4.12) to hold is constructed as 
follows. If 2h < d{u) it follows from (4.15), (4.16), and Theorem 4.1 that we may 
choose ju such that (4.3), (4.4), and (4.5) hold when V = V'^^, (j) = (j)^^, i = £{u) 
and / = f{u). 

If d{u) < 2h we simply choose 

7« = X [<t>u {-h^A - yTF) ^ 

where x is again the characteristic function of the interval (—00,0]. I.e., 7„ is the 
projection onto the non-positive spectrum of (f)u (— /i^A — V'^^) 4>u- Here we are 
considering as a multiplication operator. 
Finally, we set 

(4.23) 7 = / (t>uluMuy''du. 

Jd{u)<R 

By the properties (4.15), 7 is a density matrix. 

□ 



5. Proof of main theorem 

The proof of the main theorem on the molecular ground state energy is a rather 
standard application of the results presented in the previous Sections. We prove 
lower and upper bounds. 

Proof of Theorem 1.1. The starting point for a lower bound is the Lieb- 
Oxford inequality (2.5) from which we conclude that if '0 is a Z-particle {N = Z) 
then 

z 

{i^, H{Z, i?)^) > ^ (v., [ - i A, - V{Z, R, Xi)] V> + D{p^) - C / pf. 

i=l ■' 

In order to bound the last term we use the many-body version of the LT inequality 
(2.3). For all < £ < 1/2 we have 

1=1 ■' ■' 

Here we have used Holder's inequality for the p^/^ integral and the assumption that 
f/' is a Z-particle state. Thus 

(V,i?(z,i?)^) 

> 2Tr[ - i(l - £)A - V^^{Z, R, •)] _ - D{p^^{Z, R, ■)) - Ce-^Z. 

Here we have applied (2.8), the fact that the Coulomb kernel is positive definite 
such that D{p — p^^) > 0, and the fermionic property of the wave function. 
If we now use the scaling property (2.10) we find that 

TV[-i(l-£)A-FTF(Z,i?,.)]_ = |#/^TV[-i(l-e)|Z|-2/3A-yTF(^^r,.)]_, 
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where z = (zi, . . . , zm) and r = (ri, . . . , tm)- Using now (4.9) (with h 

and (2.18) we see that 
2Tr[-i(l-e)A-yTF(^^^^.)]_ 



= (1 - e)-^'^\Z\V^ {E'^^iz, r) + D{p^^ iz,v, •)) + (1 - s) 



-i\Z\' 



M 



124 



fe=l 



+ O(|Z|2-l/30) 

(1 - £)-3/2 (£;TF(Z, i?) + D{p^^{Z, R, .)) + i(l - e)-' Yl 



M 



fe=l 



+O(|Z|2-V30). 

We have here used the TF scaling E'^^iZ, R) = \Z\''/^E'^^ {z, r) and D{p^^ {Z, R, •) 
= \Z\'^/^D [p'^^{z,r,-)). Choosing e = completes the proof of the lower 

bound. 

The starting point for an upper bound is Lieb's variational principle, Theo- 
rem 2.4. By a simple rescaling the variational principle states that for any density 
matrix 7 on L^(]R^) with 2Tr7 < Z we have 

E{Z, R) < (2Tr [(-i|Z|-2/3A - V{z, r, x)) 7] + \Z\D{2\Z\-^p^)) . 

As for the lower bound we bring the TF-potential into play, 



\Z\-'^/^E{Z,R) < 2Tr\(-^\Z\-^/^A-V{z,r, 



,x) J 



\Z\D{2\Z\-'p^) 



= 2Tr 



[(-\\Z\-^/^A-V^^{z,v,x)) 7J 

(5.1) +\Z\D (2|Z|- - P^^(z, r, •)) - \ZW{z, r, •)). 

We now choose a density matrix 7 according to Theorem 4.2 with h = ■\/l/2\Z\~^/^. 
Note that with this choice of h we have that 

(67r2/i=^)-VTP(z,r,:c)3/2 = \Z\p^^ {z,r,x)/ 2. 

Since / p'^^{z,r,x) — J2jLi = 1 we see from (4.12) that 

2Tr7 < \Z\ + C|Z|2/3-i/i5 = 1^1(1 ^ C|Z|-i/3-i/i5). 

Thus if we define 7 = (l + C|Z|-i/3-i/i5)-i;j; .^^e see that the condition 2Tr7 < \Z\ 
is satisfied. 

Using (4.11) we conclude that 

\Z\D{2\Z\-V,-p^^{z,r,-))<C\Zf^'-'^'', 

and thus 

(5.2) |Z|D(2|Z|-V.-P^^(z,r,.)) <q^r-^/l^ 

where we have used that D (p^^(z, r, •)) < C. 

Finally, if we use (4.9), (4.10), and (2.18) we arrive at 



2Tr 



[-l\Z\-y^A-V^^{z,v, 



X) 7 



< \Z\{E'^^iz,v)+D{p^^{z,v,.)) 

+^-^f:zi+o{\z\y^-'/^^). 



fe=l 
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Since E'^^{z,r) < C and D (^p'^^{z,r, •)) < C we see that the same estimate holds 
for 7 replaced by 7. If we insert this estimate together with (5.2) into (5.1) and 
use again that E'^^{Z,R) = |Z|^/3i;TF(z,r) we arrive at the upper bound. □ 
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